Analysis of photon-atom entanglement generated by Faraday rotation in a cavity 
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Faraday rotation based on AC Stark shifts is a mechanism that can entangle the polarization 
states of photons and atoms. We study the entanglement dynamics inside an optical cavity, and 
characterize the photon-atom entanglement by using the Schmidt decomposition method. The time- 
dependence of entanglement entropy and the effective Schmidt number are examined. We show that 
the entanglement can be enhanced by the cavity, and the entanglement entropy can be controlled 
by the initial fluctuations of atoms and photons. 
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The interaction between quantized optical fields and 
atomic ensembles provides physical models for exploring 
and realizing novel applications in quantum communica- 
tion. A useful mechanism is the AC Stark shifts of off- 
resonant light with two polarizations. As the frequency 
shift in each polarization is conditioned by the num- 
ber of atoms at the corresponding magnetic sub-levels, 
any imbalance of atom numbers would lead to a rota- 
tion of Stokes parameters of light. The effect can be 
considered as a kind of Faraday rotation that correlates 
the polarization variables of photons and atoms, and it 
has been investigated theoretically and experimentally 
in the context of QND measurements and spin squeezing 
0, 0, Q ■ Interesting applications based on the interac- 
tion have also been discussed in recent literature. These 
include, for examples, quantum states teleportation and 
swapping , quantum memory [j| , generation of macro- 
scopic superposition states la, and the implementation 
of Deutsch- Jozsa algorithm [2| . 

The key to most of these intriguing applications is the 
establishment of quantum entanglement between pho- 
tons and atoms, governed by a model Hamiltonian 

Hj = 2gS z J z , (1) 

where S z is one of the Stokes operators of the field, J z is 
a collective spin operator of atoms, and g is a coupling 
constant. Such an interspecies entanglement can be fur- 
ther converted into non-classical spin states by projective 
measurement of the light field. Previous studies of the in- 
teraction (1) have emphasized the quantum correlations 
of observables relevant to QND experiments Q, 0, 0- 
From the view point of quantum information, a better 
understanding of entanglement can be gained by analyz- 
ing directly the photon-atom state vector in the Hilbert 
space. In particular, it would be useful to quantify the 
degree of entanglement and determine what parameters 
controlling the entanglement dynamics. 

In this paper we approach the problem by studying 
the Schmidt decomposition of entangled states gener- 
ated by (1). The distribution of Schmidt eigenvalues and 
the structure of Schmidt eigen-modes provide a complete 
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FIG. 1: (Color online) Atomic levels and interaction scheme 
of our model Hamiltonian. AC Stark shift is induced by non- 
resonant coupling between ground state \l) (|r)) and polariza- 
tion mode + (— ). For example, \l) and \r) refer to the two 
ground states with magnetic quantum numbers ±1/2. 



characterization of pure-state entanglement [8j- In con- 
trast to usual investigations in free space, we will examine 
the interaction (1) in an optical cavity. This is motivated 
by the fact that a strong atom-field coupling is often de- 
sirable to overcome decoherence. In free space, although 
the coupling strength can be increased at frequencies near 
resonances, the higher excitation rate would lead to sub- 
stantial decoherence from spontaneous emission 0, 0- 
The employment of a cavity provides a way to enhance 
the coupling strength while keeping the fields far away 
from resonance. In addition, the strong coupling may 
also allow us to explore cavity QED effects with fewer 
number of photons and atoms. 

Our system consists of Na four-level atoms interacting 
with a quantized cavity field with two polarizations. The 
two degenerate atomic ground states are denoted by \l) 
and |r), and they are coupled to the excited states |e;) 
and |e r ) by absorbing or emitting a photon of the corre- 
sponding polarization (Fig. 1). In our scheme, the field 
frequency is significantly detuned from the atomic transi- 
tion frequency, so that the excited states are rarely popu- 
lated. Hence the decoherence effects due to spontaneous 
decay from |ej) and |e r ) can be safely neglected. The 
off-resonant interaction leads to AC Stark shifts, which 
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are described by the effective Hamiitonian 0: 

N A 

H = gJ2 ( a + a + Mkk (l\ + \r) kk (r\) (2) 

fc=i 

with a+ and a_ being the annihilation operators of the 
+ and — polarizations of the cavity field. The coupling 
constant g = — \£Iq\ 2 / A is obtained by adiabatic elimi- 
nation of the excited levels in the large detuning limit, 
A S> |Oq|, where |Oo| is the vacuum Rabi frequency de- 
fined by the cavity mode. In our model, the atom-field 
coupling strength g is assumed to be the same for each 
atoms. A ring cavity supporting traveling wave modes 
is more suitable to realize our model. This is because 
g depends on the absolute square of mode amplitude, a 
traveling wave mode would give roughly the same value 
of g inside the ring cavity. 

By letting a ( z k) = \l) kk (l\ - \r) kk (r\, J z = ±£ fc ^*° 
and S z — ^{a + a + — a_a-), the Hamiitonian (2) can be 
expressed as: H — |iV y i(a 1 f a + + aLa_) + 2gS z J z . Since 
the first term corresponds to a constant shift of field fre- 
quency for a fixed particle number, it can be eliminated 
in a rotating frame. Faraday rotation is described by the 
second term, which is the same as (1). Noticing that 
J z and S z are formally equivalent to angular momen- 
tum operators, the corresponding eigenvectors are de- 
fined by: J z \j,m) = m\j,m), J 2 \j,m) = j(j + l)\j,m), 
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1(1 + 1)11,1). Here 
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2j = Na is the number of atoms, which is fixed in our 
problem. The n and s are difference and sum of photon 
numbers in the two polarization modes respectively. 

Initially, the system is in a product (disentangled) 
state: 



|tt(0))= | E A m \j,m)]®(J2 E P '-. 

\s—0 n- 

where A m and P Sin are amplitudes describing the initial 
state of atoms and photons respectively, and the primed 
summation sign refers to the step size equal to 2, i.e., 
n= — s, — s + 2, — s + 4, +s. Based on the Hamiitonian 
(1), the state vector at time t is given by, 
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where the phase factor e 1 9 mnt is responsible for the en- 
tanglement. Since the quantum number s plays no role 
in the phase factor, it is convenient to define a field state 
\n) f by summing s at a given n, 



Here F„ = |-P s .n| 2 is a normalization factor. With 

s=\n\ 

some careful counting of states, we note that the s sum- 
mation in Eq. (JSJ refers to the step size equal to 2 (i.e., 
s = \n\, \n\ + 2, ...). We also note that the orthogonality 
relation f(n'\n)f — <5„< jn is satisfied. 

The time-dependent state vector (4) can now be ex- 
pressed in a compact form: 

oo j 

l*W> = E E A m F n e-^ mnt \m) a \n) f , (6) 

n— — oo m—~j 

where \m) a is a short notation for the atomic state 
\j = NA/2,m). The evolution of the system is periodic 
with the period t — 2w/g. However, we shall see below 
that a significant degree of quantum entanglement can 
be established in shorter time scale, which depends on 
the fluctuations of atom (photon) numbers in the each 
atomic (polarization) states. 

To reveal the pairing structure of photon-atom entan- 
glement in the state (6), it is customary to perform the 
Schmidt decomposition so that, 



oo 

*(*)> = Y,^k\u k )a\v k )f- 
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Here X k are Schmidt eigenvalues, and \uk) a an d \vk)f 
are Schmidt eigenmodes of atoms and photons respec- 
tively. Since the Schmidt modes are orthogonal, i.e., 
a(uk>\u k ) a =j (v k '\v k )f = S kk >, if the field is found in the 
state \v k ), then with certainty the atomic states must 
be in \u k ). In addition, the distribution of Schmidt 
eigenvalues provides a measure of the degree of entan- 
glement. This can be quantified by the entanglement en- 
tropy S = — ^2 k X k In Afe, but a more convenient measure 
is the effective Schmidt number JC defined by 
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/C can be interpreted as the 'average' number Schmidt 
modes in the expansion. The higher the value of /C, the 
higher the entanglement, which shares similar features as 
the entropy. 

The Schmidt decomposition can be carried out by di- 
agonalizing the reduced density matrix of individual sub- 
systems fSj. For general state vectors, however, one needs 
to perform the diagonalization numerically. In order to 
gain analytic insight about the properties of entangle- 
ment, we consider a class of states in which A m and F n 
are defined by gaussian functions, and in addition, these 
amplitudes change smoothly with m and n: 
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Here \A and xf are normalization constants, and mo 
(no) and a a (cf) are parameters corresponding to the 
peak and width of atomic (photonic) amplitudes. In 
order to have A m described by a whole Gaussian in 
the range —Na/2 < m < Na/2, we require Na/2 > 
2a a + \mo\ so that A± Na / 2 is negligible. We point out 



that the gaussian amplitudes (9) and (10) can be em- 
ployed to capture or approximate a variety of states. 
These include atomic spin coherent states [ill and co- 
herent states of photons 0] . 

Now we make use of an identity derived from Mehler's 
formula, 



rj ' — ' 

e -* 2 /°l-y 2 /°l e -i 9 xyt = ^/i _ f(t)Y,» k Wk(x,t)V k (y,t) 

k=Q 
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where Uk(x,t) and V k {y,t) are oscillator mode functions 
that can be expressed in terms of Hermite polynomial 



Uk{x ' t} = ^ a A jk*k\ { ' l)k/2Hk {iX,CXA) 



V k (y,t) 



of ^2 k k\ 
In writing Eq. (|11|) . we have defined 
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(14) 
(15) 



with the dimensionless time r = gt/2. 

With the help of Eq. i|ll|) , we treat m — > x and n — > 
?/ and obtain an approximate Schmidt decomposition of 
the state Eq. (6). Specifically, we have the Schmidt 
eigenvalues 

A fc « (1 - "V fe (16) 

which decreases exponentially with mode index k, and 
the Schmidt eigenvectors, 

3 

\ u k) a ~ U k (m-m ,t)\m) a , (17) 

m=—j 
oo 

k)/~ I] ^(n-no.tJIn)^ (18) 

71 — ■ — OO 

where \m) a = e m °( mo " 2m ) T |m) Q and \n) f = 
e im {n -2n)T are Qjegnej^ xhe orthogonality relations 
between Schmidt eigenvectors can be satisfied approxi- 
mately: 

a(u k \ui) a = J2 U * Ul ~ jdmU* k Ui = 8 u (19) 

m=— j 

oo . 

/<«*!«!>/= E V k* V i~ JdnV k *V l = S kl , (20) 



provided that U k and T4 change smoothly with integer 
arguments. Therefore Eq. (19) and (20) are limited to 
modes with low fc's, which is sufficient at short times 
when the expansion is dominated by low modes. 

It is interesting that the widths of the gaussian factor in 
U k and V k decrease with time, due to the time-dependent 
scale factor £. Such a narrowing is not obvious because 
the left side of Eq. (11) seems to have widths that are 
constant in time. However, by taking the Fourier trans- 
form of the e~ x l a A-v /^Fg-^yt w ith respect to x, one 
can see that the transformed function is 'squeezed' in a 
variable consisting a linear combination of x and y, where 
x is the Fourier variable of x. The same feature can also 
been seen if one takes the Fourier transform with respect 
to y. As x (y) corresponds to atom (photon) number, 
its Fourier variable may be interpreted as a phase de- 
gree of freedom conjugate to the number variable. The 
narrowing found here therefore reflect a time-dependent 
numberjhase correlation developed between two subsys- 
tems [l3j. 

The entanglement entropy can now be calculated in 
closed form: 



S(t) 



In U 2 ) + In (1 - u 2 



(21) 



with the time-dependence of u given by Eq. (14). The 
effective Schmidt number has a simpler expression: 



Kit) w y/1 + (a A a FT y 



(22) 



which is an increasing function of the dimensionless time 
r. 

Eq. (|21|l and H22|) describe how the degree of entangle- 
ment increases with time. In Fig. 2, we illustrate the time 
dependence of the entanglement entropy as given by Eq. 
(|21|l and compare it with the exact results obtained from 
numerical Schmidt decomposition. We found a very good 
agreement in the time interval from zero up to a certain 
break time tb ~ l/omax, where a max = max{f7A, <tf}- 
The existence of a break time is understood from the fact 
that U k (Vfe) changes more rapidly with m (n) as time 
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FIG. 2: (Color online) An illustration of the time-dependence 
of entanglement entropy. The solid lines show the results from 
Eq. 12 1 1 , and the lines with data points are obtained by exact 
numerical Schmidt decomposition. The initial atomic state is 
an atomic spin coherent state with Na = 2<r^ atoms, and 
the fields are in coherent states. The following parameters 
are used: (a) of = 24, <ja = 3, no — and mo=0. (b) 
o F = 24, a a = 10, n = 12 and m =2. (c) o F = 24, a a = 18, 
no = and mo=0. All the curves are with the same value of 
max{<TA, of} = 24 so that break time tb ~ 0.04 is the same. 
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FIG. 3: (Color online) An illustration of time-dependence of 
effective Schmidt number. The solid lines show the results 
from Eq. 12211 . and the line with data points are obtained 
by exact numerical Schmidt decomposition. The parameters 
are: (a) of = 24, oa = 6, no — and mo=0. (b) of = 24, 
oa = 18, no = and mo=0. 

increases, and hence our continuous variable approxima- 
tion becomes invalid at long times. Indeed, at the break 
time t — tb, the width of the gaussian factor Uk or Vk 
is of the order 1, which is minimum separation between 
integers. 

In Fig. 3, we illustrate the time-dependence of the 
effective Schmidt number K(t) as given by Eq. ijiffi. We 
see that K,(t) grows almost linearly with time with the 
rate oao~f- The results also shows a good agreement with 



exact numerical calculation before the break time. We 
point out that at r = tb, K- reaches y/l + <r 2 nin ( where 
Cmin = niin{er.A, ct.f}), which is independent of the peak 
positions no and mo. 

Having discussed the entanglement in an ideal cavity, 
let us turn to systems with a leaky cavity. By inject- 
ing photons into the cavity, the question is how the out- 
put photons are entangled with atoms. According to 
the Hamiltonian (2), we notice that the interaction ef- 
fectively shift the resonant frequencies of the cavity. For 
a given atom number difference to, which corresponds to 
the eigenvalue of J z , the frequency shift is g{NA,/2 ± m) 
for the ± polarizations. We can generalize the input- 
output relations for field operators 0] to account for 
the shifts, which reads, 

a±, out (uj; m) = e j2e± ^ m >a±,in(^ m) (23) 

where 

el 20 ± ( w:m ) _ K c + iS T igm 
k c — iS ± igm 

Here we have defined 5 — uj — lu c — gN A /2, with uj c being 
the resonant frequency of the cavity, and the cavity field 
decay rate k c . Note that the detuning S has included the 
term gN A /2. 

In the bad cavity limit k c 3> gm, we find that 6± (u>; m) 
takes a simple form: 

9±(u;m)ne (u)T— ■ (25) 

provided that 3> S 2 . Here 6q(w) = tan -1 ((5/ft c ) is 
a phase angle independent of to, and the second term 
gm/K c is independent of uj. Eq. (25) indicates that if the 
input field frequency uj is near uj c +gNA/2 compared with 
the cavity width, then the output field operator would 
pick up the phase angle =F<7to/k c , in addition to 0q(uj). 
Such a result can be translated into Schrodinger picture. 
Suppose the input field is described by a single mode 
in the Fock state \N+, N-)- in and the atomic states is 
the state \m) a , then we have the transformation for the 
input-output state vectors [l5j |: 

|AT+,iV_) in |m) a -> e l29a ^ s -^ mn /^\N + ,N^ out \m) a 

(26) 

where n — N + — N_ and s — N + + N-. Apart from 
an unimportant factor e t2e o(^) s that is not responsible 
for quantum entanglement, the correlated phase factor 
e -i2gmn/K e jg p rec i se by f the same form as Eq. (6) de- 
rived in an ideal cavity. Therefore our Schmidt analysis 
in the ideal cavity can be employed by taking the interac- 
tion time t = 2/k c , which is the finite life time of cavity 
fields. In particular, if tb > 2/k c , we have 

K w y/1 + (2a A <JFg/nc) 2 (27) 

after all the input photons are fully converted into the 
output mode. 
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It is interesting to employ Eq. (27) to estimate the 
entanglement generated in free space, i.e., without the 
cavity. This may be done by replacing g with the corre- 
sponding value in free space gf, and 2/k c by the (single 
pass) interaction time tf. Apparently, g > gf due to the 
cavity enhancement would increase the entanglement for 
a fixed interaction time. A crude analysis by Kuzmich et 
al. in Ref . suggested that gf is inversely proportional 
to the volume defined by the spatial extend of the atoms 
in free space. As g is inversely proportional to the cavity 
mode volume, g/gf could be approximated by the volume 
ratio. However, a rigorous analysis in free space is diffi- 
cult because of the involvement of multi-mode dynamics 
in general 

To conclude, we perform the Schmidt analysis of the 
photon-atom entanglement generated by Faraday rota- 
tion in a cavity. We present an approximate analytic 
formula of Schmidt decomposition, which reveals the 
Schmidt mode structures and the time-dependence of en- 
tanglement. In particular we show that the initial fluc- 
tuations of and a a are key parameters to control the 
rate of change of entanglement. This work is restricted 
to the class of gaussian amplitudes defined by Eq. (9) 
and Eq. (10), which covers a range of quantum states 
of photons and atoms that have a well defined peak and 
width. In this paper we also discuss the entanglement in 
leaky systems when cavity decay rate is sufficiently large. 
The quantum dynamics beyond the bad cavity limit is 
out of the scope of this paper. This problem is difficult 
because the phase shift of output photons becomes a non- 
linear function of m when k c is comparable with g. The 
nonlinearity, which comes from the competition between 
natural cavity response function and atom-field interac- 
tion, cannot be treated by the approximation method 
described in this paper. This is an open topic for future 
investigations. 

This work is supported in part by the Research Grants 
Council of the Hong Kong Special Administrative Re- 
gion, China (Project No. 401305). 



[1] A. Kuzmich, N. P. Bigelow, and L. Mandel, Europhys. 
Lett. 42, 481 (1998); A. Kuzmich et al, Phys. Rev. A 60, 
2346 (1999); A. Kuzmich, L. Mandel, and N. P. Bigelow, 
Phys. Rev. Lett. 85, 1594 (2000). 

[2] Y. Takahashi et al, Phys. Rev. A 60, 4974 (1999); Lars 



Bojer Madsen and K. Molmer, Phys. Rev. A 70, 052324 

(2004) ; M. Takeuchi et al, Phys. Rev. Lett. 94, 023003 

(2005) ; C. Genes and P. R. Berman, Phys. Rev. A 73, 
013801 (2006). 

[3] Daniel Oblak et al, Phys. Rev. A 71, 043807 (2005). 
[4] A. Kuzmich and E. S. Polzik, Phys. Rev. Lett. 85, 5639 

(2000) ; K. Hammerer, E. S. Polzik, and J. I. Cirac, Phys. 
Rev. A 72, 052313 (2005). 

[5] C. Schori, B. Julsgaard, J. L. Sorensen, and E. S. Polzik, 
Phys. Rev. Lett. 89, 057903 (2002); T. Opatrny and J. 
Fiurasek, Phys. Rev. Lett. 95, 053602 (2005). 

[6] S. Massar and E. S. Polzik, Phys. Rev. Lett. 91, 060401 
(2003). 

[7] S. Dasgupta, A. Biswas, and G. S. Agarwal, Phys. Rev. 

A 71, 012333 (2005). 
[8] For an overview of Schmidt decomposition, see, A. Ek- 

ert and P. L. Knight, Am. J. Phys. 63, 415 (1995) and 

references therein. 
[9] K. Hammerer et al., Phys. Rev. A 70, 044304 (2004). 
[10] R. Grobe, K. Rzazewski and J. H. Eberly, J. Phys. B 27, 

L503 (1994); W.-C. Liu, J. H. Eberly, S. L. Haan and R. 

Grobe, Phys. Rev. Lett. 83, 520 (1999); R. E. Wagner, 

P. J. Peverly, Q. Su, and R. Grobe, Laser Phys. 11, 221 

(2001) . 

[11] Spin coherent states are defined by j£) = (1 + 
|£| 2 )~V J +|j, -j), see F. T. Arecchi, E. Courtens, 
R. Gilmore, and H. Thomas, Phys. Rev. A 6, 2211 
(1972). The amplitudes of states can be well approxi- 
mated by a Guassian with a a = y Na/2 and mo = 

^a(|C| 2 - 1)/2(|£| 2 + 1) when N A » 1. 
[12] For a two-mode coherent field state |a_), 



^(s/2 + n/2)\(s/2-n/2)\ 

where (s) = a+| 2 + a?_| 2 is the average total photon 
number. When (s) ^> 1, F n defined after Eq. (5) can be 
well approximated by a gaussian with <xf = 2\/ (s) and 
n = |a+| 2 - |a_| 2 . 

[13] K. W. Chan and J. H. Eberly (unpublished). 

[14] See for example, D. F. Walls and G. J. Milburn, Quantum 
Optics, Springer, Berlin (1994). 

[15] A single mode field considered here does not necessar- 
ily mean that it is purely monochromatic. In general 
a (wave packet) mode corresponds to a general super- 
position of frequency modes, so that the input Fock 
state | A/"± ) i n = fe^finlO) /y/N± \ is defined by the operator: 

6±,in = J f(ui)a± ! i n (u)dw. A similar definition can be 

written for the output mode. For the state transformation 
(26) to be valid, f(ui) has to concentrated at uj c + gNA/2 
and its width is small compared with k c . 



